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Abstract 

We propose two new sources of a substantial tensor polarization of sea partons in the 
deuteron: the diffractive nuclear shadowing which depends on the alignment of nucle- 
ons in the polarized deuteron and the nuclear excess of pions in the deuteron which is 
sensitive to the spin state of the deuteron. The corresponding tensor structure function 
b2(x,Q 2 ) rises towards small x and we predict an about one per cent tensor asymmetry 
A 2 (x, Q 2 ) = b2(x,Q 2 )/F2d(x,Q 2 ) which by almost two orders in magnitude exceeds the 
effect evaluated earlier in the impulse approximation. The both mechanisms for tensor 
structure function break the sum rule / dxbi(x,Q 2 ) = suggested by Close and Kumano. 
We comment on the impact of tensor polarization on the determination of the vector spin 
structure function gid{x, Q 2 ) for the deuteron. 



1 



The wide spread theoretical prejudice supported to a certain extent by experiment is 
that at high energies which in deep inelastic scattering (DIS) is x — > total cross sections 
cease to depend on the beam and target polarizations (hereafter x, Q 2 are the standard DIS 
variables). For instance, the SLAC E143 data || suggest, and QCD evolution analyses of 
polarized DIS support ||, the vanishing vector polarization of sea quarks in the deuteron. 
In this paper we show that sea quarks in the spin 1 deuteron treated as a proton-neutron 
bound state have a non-vanishing tensor polarization which even rises at x — > 0. [] 

If S is spin-1 operator then the spin tensor To- = ^(SiSk + SkSi — |S 2 5jfc) is the T-even 
and P-even observable. Consequently, for DIS of unpolarized leptons the photoabsorption 
cross sections <7y , and structure functions F?, (x, Q 2 ) = -^J^ — (er^ + cr^) are of the 
form = cr + cr 2 7ifc^j^fc + ••• and depend on the spin projection Sn = A onto the target 
spin quantization axis n H, |^, 0]. Hereafter we focus on the 7*-target collision axis (the 
z-axis) chosen for n, a somewhat more involved case of n normal to the collision axis will 
be considered elsewhere. The tensor spin structure function 6 2 can be defined as 

b 2 (x,Q 2 ) = ^[Ft\x,Q 2 )+Ft\x,Q 2 )-2Fi°\x,Q 2 )] 

= \ E e 2 f x[ q { f + \x, Q 2 ) + q <f\x, Q 2 ) - 2qf\x, Q 2 )] . (1) 

By the parity conservation = and F% = F 2 , = qf . Similarly, er^ and 
define the transverse, b T = 2xb\, and the longitudinal, b L = b 2 — &t, tensor structure 
functions. The precursor of the structure function b% has been discussed by Pais in 1967 0. 

Because vector polarized deuterons are always tensor polarized, in DIS of polarized lep- 
tons on the 100 % polarized deuteron the experimentally measured vector spin asymmetry 
equals Aid = 9id/F^ . Consequently, the determination of the vector spin structure func- 
tion gid only is possible if F x were known: 

gi d (x,Q 2 ) = AiM ±) (x,Q 2 ). (2) 

One must not overlook the fact that hitherto are completely unknown and all the 

deuteron data have been analyzed under the unwarranted assumption that Fid(x, Q 2 ) for the 
X A brief report on these results was presented at the Workshop on Future Physics at HERA H 
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unpolarized deuteron can be used instead of F[ ± \x 1 Q 2 ) and/or the proper average (F^) 
for the not 100% polarized target in (|2]). In order to justify this assumption and use the 
polarized deuteron data for tests of the Bjorken sum rule one needs the direct measurements, 
and the benchmark theoretical evaluations, of tensor structure functions of the deuteron. 

In the impulse approximation (IA) of Fig. la discussed ever since one must perform the 
Doppler-Fermi smearing || 



f\x,Q 2 ) = f^f w (y) 
J y 



F 2p (-,Q 2 ) + F 2n (-,Q 2 ) 

y y 



(3) 



subject to the sum rule for the smearing function 

f^f w (y) = h (4) 

« U 



where y is a fraction of the lightcone momentum of the deuteron carried by the nucleon 
(nonrelativistically y — 1 +p z /mN, where p is the nucleon's momentum in the deuteron and 
m^v is the nucleon mass). For the pure S- wave deuteron f^\y) = /«%) and b 2 (x,Q 2 ) = 0, 
the admixture of the D-wave makes f^(y) ^ f^(y) but they differ markedly only at 
values of p such that f^'(y), f^°'(y) "Cl ||- As a result, in the IA the tensor polarization 
of quarks A 2 (x, Q 2 ) = b 2 (x,Q 2 )/F 2d (x,Q 2 ) is at most several units of 10~ 4 for x < 0.5 
H |], |l0fl - Because of the sum rule (§) the tensor polarization of sea quarks vanishes at 
x — > and the widely discussed sum rule 

J dxb!(x,Q 2 ) = (5) 

holds in the IA. Capturing on the prejudice that spin effects always vanish at high energy, 
Close and Kumano [0 conjectured that tensor polarization of partons always vanishes at 
x — > and suggested that the sum rule (H) is a model independent property of DIS off spin 
1 targets. Notice, that in a disguised form the sum rule (|) first appeared in fl2|| , see also 



[I3|. For the constraints on tensor structure functions from > see higher twist 



contributions to tensor structure functions are discussed in [ p~5|j . 

In this paper, treating the deuteron as a proton-neutron bound state, we report the first 
beyond-the-IA evaluation of the tensor polarization of sea partons in the deuteron from 



Gribov's inelastic shadowing diagrams of Fig. lb and lc ||T6[ when the hadronic state X 



produced by the beam photon on one nucleon of the deuteron regenerates the photon on 
the second nucleon leading to a nuclear eclipse effect: crtot{l*d) = crtot(j*p) + crtat(j*n) — 

' s h{Y d) . There are two major polarization dependent contributions to the eclipse effect. 
The first one is the diffractive nuclear shadowing (NSH) when the state X is excited by 
pomeron exchange (Fig. lb). It is sensitive to the alignment of nucleons in the polarized 
deuteron (|| and references therein). In the second mechanism the state X is excited by 
pion exchange (Fig. lc). We comment that apart from sea quarks both mechanisms make 
gluon densities too dependent on the tensor polarization of the target deuteron (for early 
discussion on gluon densities for spin axis n normal to the collision axis see || [17], |i~8f ) 

Hereafter we focus on small x and do not consider negligible IA terms. We evaluate 
diffractive NSH using Gribov's theory |1(J] extended to DIS in ||19|| : 



7T 47T Oi em J J 



(6) 



dtdM 2 

Here M is the mass of the state X in diffractive DIS 7*iV — > XN, k = (kj_,/c z ) is the 
momentum of the pomeron, the polarization dependence comes from the deuteron form 
factor Sj^\4k 2 ). Following the conventional wisdom, we neglected the spin-flip pomeron- 
nucleon coupling. It is customary to describe diffractive DIS in terms of the pomeron 
structure function F 2 tp(P) an d the flux of pomerons in the proton 4>jp(xjp)/x-fp: 

i r>2\ da T{l* -» X) _ atotjpp) 47T 2 a em 
1 + W ' dtdM 2 ~ 16tt Q 2 

^(x!p)F 2 ^(/3,Q 2 )e-^l'l +<p s £ a (xTp)F^(f3,Q 2 )e~ B - w } , (7) 

where (3 = Q 2 /{Q 2 + M 2 ) has a meaning of the Bjorken variable for DIS on pomerons and 
xjp = x/(3 is a fraction of the nucleon's lightcone momentum taken away by the pomeron, 
in the target rest frame k z = xjpmN. We use the predictions [^TJ for structure functions 



F%j£:,F%jp and flux factors c/>^(xip), <fi^> a {x-sp) for the valence and sea of the pomeron from 
the color dipole approach to diffractive DIS |2t], |22|] which is known to agree well with the 
HI and ZEUS data [^3], |24|. In (|7|) we suppressed a small contribution from diffractive 



excitation of the open charm |20], |25| . The experimental determinations of diffraction slopes 
are as yet lacking, following |22] in our numerical estimates we take B sea = B^jp = 6 GeV~ 2 



and B va i = 2B sea , the results for NSH only weakly depend on B vai , B s 



With the input ([/]) NSH can be cast in the convolution form 



AF^\x,Q 2 ) = [ 

Jx 



1 dx 



ip 



x Xjp 



x 



X 



2 iPV— , Q 2 ) + AnW(a*)i^(— , Q 2 
Xjp 'ip 



(8) 



where An| A ^(xip) is a nuclear modification of the pomeron flux functions (i = vol, sea): 



An^ (x 



IP/ 



2 CTtotjPP) 
7T 167T 



(9) 



We recall that diffractive DIS and NSH are leading twist effects fl9| , |20|| . Consequently this 
mechanism gives a manifestly leading twist b2{x,Q 2 ). In the momentum representation 



d 3 P 



u (p')w 2 (p) 



V2 

w 2 (p)w 2 (p') f 



(2 

3(d*p)(dp)-(d*d) 



(27 



+ 



u (p)w 2 (p') 



V2 



3(d*p')(dp') - (d*d) 



9(d*p')(dp)(pp') - 3(d*p)(dp) - 3(d*p')(dp') + (d*d)] } , (10) 

where |26| x I / a(p) = (v / 2«o(p)(cd) + w 2 (p) (3(<xp)(dp) — (<rd))), cr are the Pauli spin 
matrices, p = p/|p| and p is the proton's momentum in the deuteron, p' = p + k, the S- and 
D-wave functions are related to the familiar radial functions as Uq(p) = drru(r)jo(pr) and 
w 2 {p) = J °° drrw(r)j 2 (pr), d(A = 0) = e z and d(A = ±1) = -^(T e x —i e y) are the deuteron 
polarization vectors. The alignment of the deuteron and nonvanishing contributions to 
b 2 (x,Q 2 ) come only from the S-D wave interference terms oc u (p)w 2 (p') + u (p')w 2 (p) and 
the small D-wave term oc w 2 (p)w 2 (p') . 

The pion exchange dominates 7*p — > nX, contributes substantially to 7*p — > pX at a 
moderately small x^p [B7|, and is an well established source of the u-d asymmetry in DIS 



on the nucleon [29, 30, 31, 32 . The diagram of Fig. lc for DIS on the deuteron can be 
interpreted as a contribution from DIS on nuclear excess pions in the deuteron. Here we 
focus on pions because the pion exchange is long ranged and arguably must be the dominant 
mesonic effect. The contribution of excess pions to AF s h can be cast in the convolution form 



Jx Z 



\Ff{-,Q 2 )+ l -Ff{*,Q 2 ) 
3 z 3 z 



(11) 



where F£ is the pion structure function as parameterized, for instance, in [[33|, |34| , and the 
flux function for the excess pions in the deuteron equals 



4tt 



d 2 k { Gl 



1 lI - v!J) 1 ,/>Sp Tr{*t (p^k)* w (p)(<rk))}, (12) 



(27T)2 (t-^)2 J (27T) 



where k = (kj_, k z ) is the momentum of the exchanged pion, k z = zitln (notice a similarity 
between z and xjp), t = — (k 2 L + z 2 m 2 N )/(l — z), g n NN is the uNN coupling and G^NN^^t) is 
the 7riViV form factor. Evidently, AF S ^ is of leading twist. At z <C 1 relevant to the present 
analysis the Regge model parameterization G n NN(z,t) = exp(— (Aq + a4log(-))(|t| + ji 2 ,)) 
is the appropriate one (for instance, see p5| ). With A 2 , = 0.3 GeV -2 and the slope of the 
pion trajectory o4 = 0.9 GeV -2 , the flux of pions in the nucleon evaluated with the above 
Regge form factor is for all the practical purposes indistinguishable from the result of the 
lightcone approach |29|, |30fl . In terms of the S- and D-wave functions uq,w 2 we find 



rr[^ A) (p')(o-k)* (A) (p)( CT k)] = 4tt {u (p)uo( P ') (2(kd)(kd*) - k 

u (p)w 2 (p' 



V2 
u (p')w 2 (p) 

w 2 (p)w 2 (p' 



3(d*p') (2(p'k)(kd) - k 2 (p'd)) - (2(kd)(kd*) - k 2 

3(d*p) (2(pk)(kd) - k 2 (pd)) - (2(kd)(kd*) - k 2 ^ 

[9(dp)(d*p') (2(kp)(kp') - k 2 (pp')) - 3(d*p') (2(p'k)(kd) - k 2 (p'd) 

- 3(d*p) (2(pk)(kd) - k 2 (pd)) + 2(kd)(kd*) - k 2 ] } (13) 



In the evaluation of this matrix element we differ from the expression given in [P^j . These 
authors considered deep inelastic scattering on unpolarized deuterons. Obviously, in this 
case the pure S-wave contribution must be of the form u (p)u (p')k 2 . The expression in [[36] 
does not have this form. 

Breaking of the Callan-Gross relation &l(x,Q 2 ) = is of the origin similar to that in 
the proton structure function. For instance, bi from excess pions is found replacing F^" 
in (11) by Fl(x,Q 2 ). For the evaluation of hi from diffractive NSH one needs to know 
derail* ~^ X) / dtdM 2 . The experimental data on da® are as yet lacking, here we use the 
perturbative QCD results jp"7[] , according to which R = da® I 'da® is very similar to that for 



ep and/or en scattering [B7 



All the calculations have been performed for the realistic Bonn and Paris |39| wave 
functions which are representative of the weak and strong D-wave admixture in the deuteron. 
Because of the deuteron form factor in the integrand of @ and (|T2"D An£\z), An^H^ip) 
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and An^(xp) are nonvanishing only at 

xv>,z<z(z)~-J: 0.2, (14) 

K d m N 

where is the deuteron radius. To have a quick impression on the effect of excess pions 
notice that at x < (z) the convolution ( [XT] ) can be approximated by 

The numbers of excess pions (An^) = Jq —An^'(z) are cited in Table 1. The pure S- and 
(negligibly small) D-wave terms give (An^) > 0, i.e., the antishadowing effect in AF s h- The 
tensor polarization of sea quarks is proportional to (An^) — (An^) = 1.16 • 10~ 2 , whereas 
NSH for unpolarized deuterons is proportional to | (j^An^) + (An^)^j = 0.42 ■ 10 -2 
for the Bonn wave function. Very similar results are found for the Paris wave function: 
(An^) - (An^) = 1.22 • 10~ 2 and § (2An&) + (AnjW)) = 0.34 • 10~ 2 and in Fig. 2 
and Fig. 3 we show the results only for the Bonn wave function. For a comparison, for 
the number of pions in a nucleon one finds (n n )N ~ 0.2 |29|, [30|. In Fig. 2 and Fig. 3 we 
show the contribution from excess pions to AF s h(x, Q 2 ) and 6 2 Q 2 )- All numerical results 
are for Q 2 = lOGeV 2 and for the GRV pion structure function ||33|| , we show them for the 
experimentally accessible range of x. 

The useful scale for NSH corrections at small x is set by the deuteron spin averaged 
ratios Ar^ a = An sea (xjp)/(f)^ a (x F ) and Arp = An va i(xjp)/(j)^ l (x F ), where An^Xjp) = 
I Arij- (xjp), i = val,sea . We find Arp a = 0.1, Ar^p = 0.08 for both the Bonn and 
Paris wave functions. Notice that Arj^' « Arp a because of B e i, B 3 jp -C R% which is natural. 
Our estimate of NSH for unpolarized deuterons is shown in Fig. 2. Because the pion effect 



is small, our results for NSH are close to those of Ref. |40| in which only the pure diffractive 
NSH has been evaluated in the related model for diffractive DIS. 

The results for b 2 {x, Q 2 ) f° r the deuteron are shown in Fig. 3. To a crude approximation, 
for the diffractive NSH mechanism b^(x, Q 2 ) ~ ^F s h(x, Q 2 ). The excess pions play a much 
larger role in b 2 (x,Q 2 ) than in nuclear shadowing and dominate b 2 (x,Q 2 ) at x > 0.01. At 
smaller x the diffractive NSH mechanism takes over. In Fig. 2 and Fig. 3 we also show 
our estimate for bi{x,Q 2 ). Fig. 3b shows that the tensor polarization of sea quarks in the 
deuteron A 2 (x, Q 2 ) is about 1 per cent effect and rises towards small x. 



Summary and conclusions. 

We have identified and evaluated the two new sources of a substantial tensor polarization 
of the sea distribution in deuterons. Our analysis provides a benchmark calculation of the 
tensor structure function of the deuteron at small x, very similar results must follow from 
other models provided that they give an equally good description of the experimental data 
on diffractive DIS. The results for 62(^5 Q 2 ) are stable against the admissible variations of the 
deuteron wave function. The both diffractive nuclear shadowing and pion excess mechanisms 
are of essentially nonperturbative origin but give the leading twist tensor structure function 
b2{x, Q 2 ). To this end, the pion contribution is obviously GLDAP evolving, the sea structure 



function of the pomeron also obeys the GLDAP evolution \T2\, the detailed form of scaling 
violations in the valence structure function of the pomeron at (3 — > 1 remains an open issue 
p5| , |4l|j but the contribution from this region of (3 to b2{x,Q 2 ) is marginal. Consequently, 
QCD evolution of b2(x,Q 2 ) must not be any different from that of ^(x, Q 2 ), see also the 
related discussion in [|, [0| . 

Our finding of the rising tensor polarization of sea quarks in the deuteron invalidates the 
assumptions invoked in the derivation of the Close-Kumano sum rule (|5|) and we conclude 
that the sum rule (||) is faulty. 

The present uncertainties in the wave function of the deuteron, structure function of the 
pion and the diffractive structure function of the proton do not affect our principal conclusion 
of a substantial, ~ 1%, tensor polarization of sea quarks in the deuteron, which is within 
the reach of future experiments with polarized deuterons. Regarding the measurements 
of gid in polarized ed } fid scattering, our results justify the substitution of the unpolarized 
structure function of the deuteron for the more correct but still unknown F 2 in (@) at the 
present accuracy of measurements of the spin asymmetry A^, but corrections for the tensor 
polarization effects will eventually become important for the quantitative interpretation of 
high precision DIS off polarized deuterons in terms of the neutron spin structure function 
gi n and for high accuracy tests of the Bjorken sum rule. 

Finally, eqs. (|8]) and (|13|) are fully applicable to calculation of the tensor polarization 
of gluon densities in the deuteron polarized along the 7*6? collision axis in terms of gluon 



S 



structure functions of the pomeron and pion. The results for the gluon structure functions 
will be presented elsewhere. 
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figure captions 

• figure 1: Impulse approximation (a) and Gribov's inelastic shadowing (b,c) diagrams 
for DIS on the deuteron. 

• figure 2: (a) Nuclear shadowing correction in the deuteron AF sh (x,Q 2 ) for Q 2 = 
lOGeV 2 and its diffractive NSH and pion excess components for the Bonn wave func- 
tion of the deuteron. (b) The fractional shadowing correction to the deuteron structure 
function F 2d (x, Q 2 ) for Q 2 = lOGeV 2 . 

• figure 3: (a) The tensor structure function of the deuteron b 2 (x, Q 2 ) at Q 2 = lOGeV 2 
and its diffractive NSH and pion excess components for the Bonn wave function of 
the deuteron. (b) The tensor polarization A2(x,Q 2 ) of quarks and antiquarks in the 
deuteron at Q 2 = lOGeV 2 . Also shown (by the dotted line) is the longitudinal tensor 
structure function bi(x,Q 2 ) (in (a)), and the corresponding asymmetry (in (b)). 
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(AnW) • 10 2 


Bonn Paris 




A = 


A = ± 


A = 


A = ± 


s-wave 


0.146 


0.999 


0.133 


0.929 


s/d-int. 


-0.537 


-0.268 


-0.676 


-0.331 


d-wave 


0.035 


0.085 


0.070 


0.153 


total 


-0.356 


0.807 


-0.466 


0.750 



Table 1: The S-D wave decomposition of the mean multiplicities of excess pions in the 
deuteron for different polarizations A along the 7* ' d-collision axis. The numerical results 
were obtained for the Bonn and Paris wave functions. 
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